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ABSTRACT 

Radio interferometers can measure the full polarization state of incoming 
waves by cross-correlating all combinations of two orthogonal polarizations at 
each antenna. The independent sets of electronics used to detect the two po- 
larization states will introduce a differential instrumental delay between the two 
data streams. The usual technique of separate calibration of the parallel-hand 
sets of visibilities still allows for an arbitrary offset in group delay and phase 
between the two parallel systems. In order to use the cross-polarized visibilities, 
this instrumental offset must be determined and removed. This paper describes 
one such technique and explores its application in the Obit package. The tech- 
nique is successfully applied to some EVLA data using both strongly and weakly 
polarized calibrators. 

Subject headings: Data Analysis and Techniques 



1. Introduction 

Radio interferometry provides a powerful tool for probing the structure of the polarised 
emission from celestial sources. Such measurements are made by detecting two orthogonal 
polarization states at each antenna and cross-correlating all combinations. To utilize such 
observations, corrections must be applied for the effects of the independent electronics in 
addition to paths through the atmosphere and errors in the assumed geometric model that 
corrupt the data. Most of the geometric and atmospheric effects modifying the arrival time 
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of a given wavefront are the same for the two polarization states at a given antenna but dif- 
ferences in the electronics used for the two polarization states will introduce an instrumental 
polarized delay and phase. 

The traditional calibration technique for radio interferometric data, especially when the 
receptors measure the orthogonal circular polarizations, has been to determine the calibration 
for the two hands of polarization separately. This effectively ignores both the polarization 
of the calibrator and the spurious polarized response of the instrument. 

For systems measuring circular polarization and using synchrotron sources for calibra- 
tion, this is generally a good approximation. The circular polarization of the typical compact 
radio source is less than 0.1% and the spurious instrumental polarization response of the or- 
der of a few This approximation is less valid for interferometers measuring orthogonal 
linearly polarized signals as the same calibrators typically have linear polarizations of a per- 
cent or less a t frequencies below a few hundred MHz up to ~10% at frequencies above a few 



tens of GHz (ITucci & Toffolatti 1 120121 ) . 



Independent calibration of the parallel-hands also has the disadvantage that it allows 
an arbitrary group delay and phase offset between the two parallel systems. The offsets 
between the parallel-hands need to be determined and removed before the cross-polarized 
visibilities can be used for astronomical imaging. These cross-polarized visibilities may be 
dominated by the linear polarization of the source which can be used for cross-polarized 
calibration; this case is considered in the following. 

One technique used in the past was to "frin ge fit", fit for group dela y and delay rate, 



independently for the parallel-hand visibilities (ISchwab fc Cotton I Il983l ). apply this cab 
bration and then do sin gle baseline d elay fitting of a few cross -polarized visibility spectra 
and average the results J Cotton Ill99i ). An alternate technique (IBrown. Roberts fc Wardle 



19891 ) is to use closure constraints on the results of a series of single baseline fringe fits to 



determine the offsets between the parallel hand systems. Both of these techniques have lim- 



ited sensitivity. This paper explores an improved technique in the Obit package ((Cotton 
20(3)1. 



1 This is discussion is limited to tracking dishes, aperture arrays such as arrays of dipoles used at low 
frequencies can have very large and variable instrumental linear polarizations 

2 http://www. cv.nrao.edu/~bcotton/Obit. html 
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2. Interferometric Polarimetry 

Telescopes using heterodyne electronics as are commonly used in interferometers at radio 
wavelengths are sensitive to a single polarization state of the incoming wave. In order to fully 
sample both polarization states of the celestial signals, sets of electronics nominally sensitive 
to orthogonal polarizations are used. To sample the full state of the visibility measured with 
an interferometer, cross-correlations of all (4) combinations of the polarizations states are 
made on each baseline, or pair of antennas. Denote the various visibilities measured between 
antennas j and k and using detectors for polarization states p and q as v^. Typically, either 
right- and left-hand circular polarizations ("R", "L") or orthogonal linear polarizations 
( "X" , "Y" ) are used in radio interoferometers. 

In practice, the two polarizations measured are not precisely those desired but can be 
modeled by the desired state plus a complex value, called the "leakage" term, dj P , times 
the orthogonal state. Thus, the signal received by a detector on antenna j nominally in 
polarization state p is actually: 

S jp = S jP + djpSjq (1) 

The leakage terms will result in a spurious polarized response which is referred to as "in- 
strumental polarization" in the following. 

To first order (ignoring products of source and instrumental polarization terms), the 
interferometric response for an interferometer between antennas j and k using linear detectors 



of an unresolved source is ( iThompson. Moran fc Swenson 1120011 ): 
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where gj p is the complex gain of the electronics for polarization p on antenna j, "*" denotes 
the complex conjugate, /, Q, U, and V are the Stokes parameters of the source emission, % 
is \J— 1 and, x is the parallactic angle given by: 

i / cos A sin h \ 

* = tal1 — 1 X \~- A h ) 3 

Vsm A cos o — cos A sin o cos h) 



where 5 is the source declination, A is the latitude of the antenna and h is the hour angle of 
the source. For linearly polarized detectors with detectors rotated from the local horizontal 
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and vertical, this rotation needs to be added to the value of x given above. To simplify the 
notation in the following, we assume nearly identical parallactic angles at all antennas at a 
given time, although the technique is not limited to this case. 

The corresponding first order relations for antennas with circular detectors is: 

1 

vrr = ^9jRg k R{I + V), 

vrl =lg j R9*k L [(d j R + dl L )I + e~ 2lx (Q + iU)}, 



(4) 



T 

vlr =\g ] Lgl R [{d 3L + dl R )I + e 2ix (Q - iU)), 

vll = ^gjigldi - V). 

The impact of source polarization on parallel hand calibration can be seen from eqs (2) and 
(4). Stokes V is generally very small compared to Stokes I, much less than a percent and 
can relatively safely be ignored in eq (4). On the other hand, values of Stokes Q and U of 
a few percent of Stokes I will limit the accuracy of parallel hand calibration if ignored in eq 
(2). 



Sault. Hamaker fc Bregman I (119961 ) give a discussion of the three calibration parameters 
that cannot be constrained from only parallel hand data. If the source and instrumental 
polarizations can be considered negligible, the third and its frequency derivative can be 
obtained using the linear polarization of the calibra tor. For more details on the response of 



an interferometer to a partially polarized signal see iThompson. Moran fc Swenson I (120011 ). 



3. Independent Parallel Hand Calibration 

If the calibrator and instrumental polarizations are sufficiently small to be ignored, the 
two parallel polarized sets of visibilities can be assumed to be independent measurements of 
the same celestial quantities. For a given baseline, j — k, the measured correlations {v pp ^ hs ) 
are related to the calibrated visibilities (v^l' cal ) by: 

pp-cal PP-obs p p* /r\ 

jk — u jk 9j gk \°) 



where g? is given by: 



p p -2-Ki(&T V v)-iS p 

cf\ = cij e 



and a? is the amplitude correction, ArJ the group delay residual from the correlator model, 
<p p - the model phase residual, and v the observing frequency. Calibration parameters a p -, ArJ 
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and <f) P j can be determined from observations of (calibrator) sources of known brightness, 
structure and position. Calibration quantities are generally functions of time and frequency. 

This system of equations is degenerate in that only differences in ArJ and <j) p are actually 
measured. The system is frequently made determinate by assigning the values of a "reference 
antenna" to zero. Thus, all phase-like quantities are determined relative to the reference 
antenna leaving only the relationship between the two hands of polarization at the reference 
antenna undetermined. If the electronics of the reference antenna are sufficiently stable, the 
differences between the two parallel hands at the reference antenna will be stable in time. 



4. Residual Cross— hand Group Delay/Phase Offset 

Independent calibration of the parallel-hand systems allows for a difference between 
At? and Ar q as well as between c/% and where r indicates the reference antenna. Denote 
these differences as Ar^ 9 and (j^, q . Note: if the observing bandwidth is divided into multiple 
"spectral windows" with at least partially independent electronics and signal paths, the cal- 
ibration parameters may vary between spectral windows. The relation per spectral window 
between parallel-hand calibrated data {v p £~ m ) and cross-hand calibrated data (v?l~ xca ) is: 

v pq^ccal _ v pq.cal e -2mAT? q v-i<t>P q /g\ 
jk jk V / 

Note: this requires a single pair of values, Ar^ q and <$. q , per spectral window. A sample 
baseline of EVLA (circular polarization) calibrator data with parallel-hand but not cross- 
hand calibration applied is shown in Figure [TJ 



5. A Method of Determining the Cross— hand Offset 

While there are but two parameters to be determined per spectral window, there are 
a number of complications. The most serious of these is spurious instrumental polarization 
(dj P in Eqs. |2]and|4j) which can contribute a significant fraction of the strength of the source 
polarization to the cross-polarized visibilities. To first order, this spurious instrumental 
polarization is independent of source polarization and can vary strongly with frequency and 
baseline. 

A second complication is for arrays with alt-az mounts for which the antenna rotates with 
parallactic angle as seem by the source (Eq. [3]). If no correction is made and the antennas 
detect circular polarization, a varying parallactic angle will cause the source polarization to 
rotate in polarization angle while the instrumental polarization is constant. If the data are 
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corrected for parallactic angle, the instrumental polarization then rotates while the source 
polarization is constant. 

A further complication is the case in which the source has significantly resolved polarized 
structure. In this case, the source component of the cross-polarized response will vary with 
time, frequency and baseline. 



5.1. Circularly Polarized Detectorss 

Consider the case of the EVLA with alt-az mounts and detectors sensitive to circular 
polarization observing an unresolved source whose linear polarization is much stronger than 
the instrumental polarization. If these data have their phases corrected for parallactic angle 
and have parallel-hand calibration applied, the cross-polarized spectra are from Eq. HJ 

v RL.cal ^ fg + i jj\ e +2TTi(AT? q u)+itf q 

y LR-cal ~ (Q _ i jj^ e -2-ni{AT? q u)-i<t > p r q ^ > 

Thus, the RL and LR spectra are the complex conjugates of each other and this relationship 
holds for all baselines and times. In the case represented here, the cross-polarized spectra 
can be averaged over all baselines and times and the averaged spectra used to fit for At*" 3, 
and (jf r q . 

The solution for Ar*" ? and can be separated and done independently. The group 
delay (Ar^ 9 ) can be determined from the value of At which maximizes: 

\Z{vf k L - Cal + v LR-cal*y-2ni(ATv) | ( g ) 

using a direct parameter search over At and 1 1 denotes the modulus. The summation is over 
spectra for all included times and baselines. 

The phase (4>^. q ) is then the phase of the average of the delay corrected visibilities: 

= arg[E{vf k L - ml + vf k R - cal *)e- 2wi{Ar " qu) ] (9) 

The analysis above is only valid if the cross-polarised visibilities include a significant de- 
tection of the polarized signal. The typical "signal-to-noise ratio" (SNR) of the individual 
channel samples can be estimated from the RMS scatter of their phases about the mean in 
radians and the approximation that: 

RMS ph 
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for an SNR > 1. Channel averages over baseline and time can be used to increase the SNR 
of individual samples. 

In general, the instrumental polarization is not negligible compared to the source polar- 
ization but should vary from baseline-to-baseline. Furthermore, after the data are corrected 
for the parallactic angle, the instrumental polarization will vary with parallactic angle, hence 
time. Averaging over baselines and time will reduce the contribution of the instrumental 
polarization to the averaged RL and LR spectra. 



5.2. Linearly Polarized Detectors 

The case of interferometers, such as ALMA or the ATCA, whose elements are sensitive 
to linear polarization, is more complex. As can be see from Eq. [2j the source polarization 
component of the cross-hand visibility spectrum is: 

v XY.cai K ^Q s i n2x + Ucos2x + iV)e +2ni{AT " qu)+i ^ q 

v YX.cai ^ (_Q s in2x + f/cos2x - iyy-^W 1 ^)-^ . ( ' ' ' 

Stokes V for most sources is very small and can be ignored. However, the linearly polarized 
component is a more complex function of parallactic angle and can only be averaged over a 
time range and set of baselines for which the variation of parallactic angle is small. 

The averaged XY and YX spectra can be used to determine the X-Y delay difference as 
was done for circular feeds in the section above and the SNR of the data can be determined 
as well. Variation of the instrumental polarization among baselines will reduce its effect on 
the averaged spectra. 



6. RLDly: An Obit Implementation 



The technique outlined above has been implemented in the Obit flCotton 1120081 ) package 
available through task RLDly. The fitting is that described in the previous sections except 
using weighted sums to fit for At? 9 and § v T q . The weights used are the data quality weights 
assigned to the data which are inversely proportional to the visibility noise variance. The 
implementation in RLDly allows for dropping the end channels in each spectral window which 
may be poor representations of the values being fitted. In this case the phase ramp applied 
to the spectra still begins with the first channel in each spectrum. 
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7. Testing 

The calibration data from an EVLA observation was used to test this technique. The 
data included the bandpass 6.0-8.0 GHz in 16 spectral windows (IFs) each averaged to 16 
channels. The observations included 5x15 second scans over roughly 5 hours on the strongly 
linearly polarized calibrator 3C286 (~ 10 %). The EVLA has detectors sensitive to circular 
polarization and alt-az antenna mounts. 

Data had strong, non-celestial signals removed and were initially calibrated using the 
standard parallel-hand techniques; a sample of these data on one of the longer baselines is 
shown in Figure [TJ Corrections were applied for parallactic angle so all times were used in 
the fitting. 

The data were obtained in the "A" configuration for which the synthesized beam is about 
0.25" FWHM. 3C286 is well resolved in this configuration but the bulk of the emission, both 
total and polarized intensity is from the marginally resolved core. Therefore, all baselines 
were used for the fitting. A target R-L phase of 66° and rotation measure of rad/m were 
provided. The outer 3 channels on each end of each spectral window were omitted from the 
fitting of the 3C286 data. 

The cross-polarized corrections were applied to the data. The same data shown in Figure 
[1] after this recalibration are shown in Figure |2j The rapid variation of the cross-polarized 
phase with frequency is largely eliminated. Variations of phase with frequency remaining 
are due to the frequency variable instrumental polarization which has not been corrected; 
these can be corrected using AIPS task PCAL. The fluctuations of instrumental polarization 
are also visible in the cross-polarized amplitudes. Note that the parallel-polarized data are 
unaffected by the cross-polarized calibration. 

In order to test the effect of a weakly polarized calibrator, the same test was run but 
using as calibrator J1504+1029 whose polarization (« 1 %) is less than the typical baseline 
instrumental polarization. A plot of the data shown in Figures [T] and |5] after calibration 
using this source is given in Figure |3j 

8. Discussion 

A relatively general and robust method for fitting instrumental cross-polarized group 
delays and phases is presented and tested using EVLA data. In these tests, the large phase 
slope in the cross-polarized spectra due to the difference in the parallel-hand group delays 
is essentially removed. 
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The fitting of data from circularly polarized detectors and a strongly linearly polar- 
ized unresolved calibrator, as in the first test presented, is relatively straightforward as all 
baselines are measuring essentially the same value and much averaging can be incorporated. 
More extensive testing will be needed in the case of weakly polarized calibrators where the 
cross-polarized visibilities will be dominated by instrumental polarization and arrays using 
linearly polarized detectors for which time averaging need be more limited. One test on the 
same dataset as discussed above using a calibrator whose polarized emission is less than the 
typical baseline instrumental polarization resulted in cross-polarized delay fits comparable 
to those derived from 3C286 although with larger scatter in the R-L phase. 

Fitting for cross-polarized group delays is needed prior to fitting for instrumental po- 
larization if any averaging in frequency is to be used in the calibration fitting. Fitting for 
cross-polarized phases at this point will be only approximate due to corruption by the in- 
strumental polarization and will need to be redetermined after removal of the instrumental 
polarization. 

The technique described here can also be applied to VLBI data for which the parallactic 
angle is different amoung the various antennas at a given time if the data are initially cor- 
rected for parallactic angle (circular feeds are generally used for VLBI). At VLBI resolutions 
there is generally significant resolution of the polarized emission of even calibrators and a 
limited range of baselines may need to be used. 

The author would like to thank the anonymous reviewer for many helpful comments, 
corrections and suggestions that led to a significant improvement in the paper. 

Facilities: VLA. 
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Fig. 1. — Sample EVLA visibility spectra for 3C286 averaged over a 15 second scan, parallel- 
hand visibities above, cross-hand below. Parallel-hand calibration has been applied but not 
the differences between the R and L systems. The upper plot in each panel gives the phase 
in degrees and the lower, the amplitude in Jy or mJy Note the rapid variation of cross- 
polarized phase with frequency indicating a large delay offset. 
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Fig. 2.— Data shown in Figure [D but with the R L delay and phase offsets applied. 
The strongly polarized source 3C286 (data shown) was used for the calibration. Note: 
the parallel-hand data is unaffected. Instrumental polarization corrections have not been 
applied. 
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Fig. 3. — Data shown in Figure |2] but with the R L delay determined from the weakly 
polarized calibrator J1504+1029. 



